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ABSTRACT 

This paper reports the results of a survey of Doppler shift oscillations measured during solar flares in emission 
lines of S XV and Ca XIX with the Bragg Crystal Spectrometer (BCS) on Yohkoh. Data from 20 flares that 
show oscillatory behavior in the measured Doppler shifts have been fitted to determine the properties of the 
oscillations. Results from both BCS channels show average oscillation periods of 5.5 ±2.7 minutes, decay 
times of 5.0±2.5 minutes, amplitudes of 17.1 ± 17.0 km s , and inferred displacements of 1070± 1710 
km, where the listed errors are the standard deviations of the sample means. For some of the flares, intensity 
fluctuations are also observed. These lag the Doppler shift oscillations by 1/4 period, strongly suggesting that 
the oscillations are standing slow mode waves. The relationship between the oscillation period and the decay 
time is consistent with conductive damping of the oscillations. 

Subject headings: Sun: corona — Sun: flares — Sun: oscillations — Sun: X-rays, gamma rays 



1. INTRODUCTION 

The solar corona is a high-temperature, fully-ionized 
plasma structured by magnetic forces. This structuring results 
in the familiar loop-dominated appearance of the corona. The 
fact that the corona is threaded by a magnetic field also means 
that it is subject to a rich assortment of oscillatory modes, 
which have b een the subject of consider a ble th eoretical inves- 
tigation (e.g..lRoberts et al.ll983lU984T) . lRobertsl(l2000h . and 
iRoberts &*N akariakov (2003) provide recent reviews of the 
theoretical state of the subject. 

Searches for coronal oscillations have a long history, with 
much of the earliest work centered on looking for evidence 
of the propagation of the wav es that were tho ught to heat the 
corona (for a review see, e.g., Aschwanden 2003). But it has 
been the more recent observations from spacecraft imaging 
instruments, especially the Transition Region and Coronal 
Explorer (TRACE) that have begun to expose the rich vari- 
ety of oscillatory phenomena present in the corona. TRACE 
observations of spatial shifts of 1 MK loops in a flaring ac- 
tive region showed unambigu ously that periodi c oscillations 
were occurring, which Aschwanden et al. ( 1999) identified as 
MHD kink mode oscillations. 

Spectroscopic instruments have also detected coronal os- 
cillations. Observations obtained with the Solar Ultravio- 
let Measurements of Emitted Radiation (SUMER) spectrom- 
eter on SOHO, show the presence of damped oscillatory 
Doppler shifts in the emission lines of Fe XIX at 1118.1 A 
and Fe XXI at 1354.1 A (TWang et al.l2002HKliem et al.l2002l) . 
These lines are formed at a bout 8 and 10 MK, respec- 
tively. Qfman & Wang l (|2002) interpreted these observations 
as evidence for rapidly-damped slow-mode standing acoustic 
waves. Additio nal evidence for this conclusion has been pro- 
vided by Wang et al. (2003b) and an extensive com pilation of 
SUMER observations is contained in Wang et al. (2003a). 

Both the oscillations observed with TRACE and those 
observed with SUMER generally damp in a small num- 
ber of periods, and ther e has been considerab l e spec- 
ulation about the cause {N akariakov et a ll 11999^ [ Qfman 
l2^02tlOfman & Aschwandenl2002HSchriiver & Brownl2000t 
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Ofman & Wang 2002). Additional observations, particularly 
if they cover different oscillation periods or different solar 
plasma conditions than those already measured may provide 
further constraints on the damping mechanism, possibly lead- 
ing to new insights into how the corona is heated. 

Detection of Doppler shift oscillations with SUMER in 
emission lines formed at flare temperatures suggests that other 
instruments designed to observe Doppler shifts in flaring plas- 
mas may be able to det ect and characterize these events. 
Recently, Mariska (2005) reported the detection of damped 
Doppler shift oscillations in spectra obtained with the Bragg 
Crystal Spectrometer (BCS) on Yohkoh. In this paper, I report 
on a study of Doppler shift oscillations in a large sample of 
the flares observed with that instrument. 

2. BCS OBSERVATIONS 

The BCS flew on the Japanese Yohkoh satellite, taking use- 
ful data from 1991 October 1 to 2001 December 14. Four bent 
crystals viewed the entire Sun over narrow wavelength ranges 
centered on emission lines of Fe XXVI, Fe XXV, Ca XIX, and 
S XV, with the two Fe bands sharing one detector and the 
Ca and S bands sharing a second. Useful data were rarely 
recorded in the Fe XXVI band, and the spectrum observed 
in the Fe XXV band is complex — consisting of contributions 
from several stages of ionization of Fe. Thus, this study 
was restricted to data from the Ca XIX and S XV wavelength 
bands. Furt her details on t he cha racteristics of the BCS are 
provided by Culhan e et alJ (119911) . 

To search for Doppler shift oscillations, plots of the total 
count rate in the BCS Ca XIX channel as a function of time 
for the first year of operations, 1991 October 1 through 1992 
September 30, were scanned for candidate flares. Any flare 
with peak count rate in the Ca XIX channel of more than about 
1000 counts s~ 1 and a smooth profile suggesting that the event 
was not seriously contaminated by emission from other events 
on the solar disk was processed with the standard BCS fit- 
ting software and plots of the derived physical parameters as 
a function of time were examined visually for evidence of os- 
cillatory behavior. 

In addition to flares observed during the first year of 
Yohkoh observations, all the flares examined in the study by 
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FIG. 1. — BCS spectra in the S XV and Ca XIX channels obtained on 
1991 October 21 at 2:48 UT. Fluxes are in units of photons cm -2 s A -1 . 
The data are plotted as histograms and the best-fit model as a smooth line. 
Note that the wavelength axis for both plots is based on the nominal Yohkoh 
pointing and assumes the flare is near the equator. Thus the absolute values 
of wavelength are not precise. 



iMariska & McTiernanl ([1999) of occulted and nonocculted 
limb flares were included. TRACE began routine observations 
in 1998 May, and so had some ove rlap in time with the Yo hkoh 
mission. All the flares l isted in ISchriiver et al.l (120021) and 
lAschwanden et all |2002) that were also observed with the 
BCS were therefore also included in this study. The Yohkoh 
mission also overlapped in time with the SOHO mission. A 
search of the BCS data, however, show s no useful observ a- 
tions for any of the events observed by Wang et al. ( 2003a). 

Figure^shows examples of typical spectra obtained in the 
S XV and Ca XIX channels for a flare observed at roughly 
02:48 UT on 1991 October 21 along with theoretical spec- 
tra fitted to the data. Although only a few strong features 
are obvious in both channels, the spectra are actually quite 
complex. Thus the fitting is done by using detailed atomic 
physics data for the transitions that produce lines in the 
wavelength ranges of the two channels. A t omic dat a for 
the fits are primarily from iBelv-Dubau et all lll982blal) and 
IVainshtein & Safronoval (H978T 1985). The required ioniza- 
tion eq uilibrium calculations are from lArnaud & Rothenflud 
(1985). 

The S XV and Ca XIX resonance lines, which are the 
strongest lines in each channel, are formed over relatively 
broad temperature ranges. The contribution functions peak 
at about 15.8 and 31.6 MK, respectively. In typical flares, 
however, most of the emission comes from plasma with tem- 
peratures less than about 20 MK, so the emission observed 
in the BCS Ca XIX band tends to be from plasma at lower 
temperatures than the peak of the contribution function. The 
key feature of the emission in both bands, is that there are 
line ratios in each band that are temperature sensitive. Thus, 
carefully fitting all the emission lines in each band provides 
an excellent measurement of the temperature of the emitting 
plasma. 



Each spectrum is fitted using an isothermal model and the 
standard BCS fitting software. Thus the fit is characterized 
by a temperature T, an emission measure EM, a nonthermal 
broadening and a Doppler shift velocity v. The best-fit 
values for the parameters are determined using Levenberg- 
Marquardt least-squares minimization (e.g., Bevington 1969). 
For the S XV spectrum shown in the figure, the best fit param- 
eters for T, EM, and % are 12.6 MK, 2.30 x 10 48 cm" 3 , and 
38.1 km s _1 , respectively. For the Ca XIX spectrum, the cor- 
responding values are 13.6 MK, 2.30 x 10 cm~ 3 , and 76.7 
km s _1 . As Figure ^shows^ the fits are generally excellent. 
The different values for T and | suggest, however, that the 
two wavelength channels are not sampling exactly the same 
plasma. 

The BCS is uncollimated and therefore views the entire 
Sun. When only one flare is taking place, this does not present 
any problems, since the high-temperature plasma is generally 
concentrated in a small enough volume that very little line 
broadening due to the spatial extent of the plasma is present. 
Because the angle of incidence of the emission into the spec- 
trometer depends on the flare's location, the reference wave- 
length for each flare will be different. The wavelength scale 
shown on the spectra in Figure ^ is f° r the nominal space- 
craft pointing and a flare near the equator. Since this study 
is concerned with Doppler shift variations within each indi- 
vidual flare, the lack of an absolute reference wavelength is 
unimportant. 

In its normal mode of operation, the BCS obtains spectra in 
all four wavelength channels every 3 s. Early or late in flares, 
when the count rates can be quite low, the individual spectra 
are difficult to fit. For this work, I have therefore summed 
individual spectra until a minimum of 10,000 counts are ac- 
cumulated in each channel. This generally results in excellent 
fits to the data and assures that each data point in the resulting 
time series will have roughly the same statistical significance. 
Note that these accumulated spectra are produced for each 
channel individually. Thus the time series of fitted spectra 
for the two channels will not generally be at identical times. 

BCS data for a total of 103 flares were processed in the 
manner described above and examined for evidence of oscilla- 
tory behavior. Of this total, 38 showed evidence of oscillatory 
behavior in the measured Doppler shifts that looked promis- 
ing enough to analyze further. Of those 38 flares, only 20 
were deemed suitable for inclusion in this paper. While the 
remaining 18 flares could be fitted using the functional form 
developed in the next section, the resulting fits were less con- 
vincing than the 20 included flares. Generally, this was due to 
the Doppler-shift data exhibiting less than a complete oscil- 
lation period, leading to considerable ambiguity in the decay 
time determination. 

3. RESULTS 

Figure|5]shows the temporal behavior of the count rates and 
the fitted parameters from the fits to the accumulated spectra 
for the S XV and Ca XIX wavelength channels for one of the 
flares that exhibited damped oscillatory behavior. The count 
rates, temperatures, emission measures, and nonthermal ve- 
locities show relatively smooth behavior. In both channels, 
however, the Doppler shifts show clear evidence for oscilla- 
tory behavior. In particular, the Doppler shift as a function of 
time in the S XV channel appears to show clearly evidence for 
a damped oscillation. 

Both the BCS S XV and Ca XIX wavelength channels are 
within the GOES 1-8 A band, and the S XV intensity curve 
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FIG. 2. — Temporal behavior of the intensity, temperature, emission measure, nonthermal broadening velocity, and Doppler shift derived from the BCS S XV 
and Ca XIX observations for the 1991 October 21 flare. The zero value for the Doppler shift velocities has been set to the wavelength shift of the first fitted 
spectrum. 



shown at the top left of Figure [2] generally tracks the inten- 
sity measured with GOES in that band. For example, the peak 
flux in the GOES 1-8 A band occurred at 2:50 UT, the same 
time as the peak in the BCS S XV band. The emission in the 
BCS Ca XIX wavelength channel tends to more closely track 
the GOES 0.5^4 A band behavior. Thus, the emission in the 
BCS Ca XIX wavelength channel peaks earlier than that seen 
in the S XV channel. As both the S XV and Ca XIX inten- 
sity and Doppler-shift data show, oscillatory behavior is gen- 
erally present at or before the peak intensity of the flare and 
is no longer detectable later in the decay phase. For flares in 
which the entire rise-phase is captured, the oscillations are not 
initially present, but do develop as the count rates approach 
maximum. 

Since the BCS is uncollimated, the location on the Sun of 
each flare determines where on the detectors the spectrum 
falls and there is no absolute wavelength reference. In the 
past, most analyses of Bragg crystal spectrometer data have 
adopted the practice of setting the rest wavelength by as- 
suming that late in the flare Doppler shifts have disappeared. 
Since a cooling plasma can often result in downflows, this 
approach is not without possible ambiguity. I have therefore 
chosen to set the rest wavelength using the first spectrum in 
each set of flare observations. Examination of the velocity 
plots for the BCS S XV and Ca XIX wavelength channels in 
Figure |2] shows that for this flare, assuming that the zero for 
the velocity scale can be set using a measurement late in the 
flare, suggests that early in the event the measured velocities 
represent roughly 10 km s _1 upflows. The relationship be- 
tween these small Doppler shifts and the much larger Doppler 
shifted secondary component that appears very early in many 
flares is not clear. Thus considerable caution should be exer- 
cised in drawing any connection between the overall trend in 
the Doppler shifts and a possible driver for the oscillations. 



Following the example of the ea rlier work on spatial oscil - 
lations observer with TRACE (e.g. JAschwanden et alJ l2002). 
I fit the BCS Doppler shift observations with a combination of 
a damped sine wave and a polynomial background. Thus, for 
each channel, I assume that the data can be fit with a function 
of the form 

v(t) =A o sin(a>r + 0)exp(-Ar)+fi(f), (1) 

where 

B(t) =b + bit + b 2 t 2 + b 3 t 3 + ■■■ (2) 

is the trend in the background data. Beginning with an initial 
guess for the fit parameters, the data in each channel were fit- 
ted to equation Q u sing Leve nberg-Marquardt least-squares 
minimization (e.g., Bevington 1969). For all the flares stud- 
ied the number of terms needed for the background expression 
was generally one or two, although occasionally three terms 
produced a better fit. 

While the S XV Doppler shift data were generally quite 
smooth, often the Ca XIX data showed considerable fluctua- 
tions on top of the apparent damped oscillation. These fluctu- 
ations are present even though the data have been accumulated 
to more than 10,000 total counts in the channel, and are proba- 
bly a manifestation of the fact that not all the high-temperature 
flare plasma is taking part in the oscillation. In some cases 
these secondary fluctuations were sufficiently large that it was 
difficult to fit the time series using equation Q. For the flares 
that showed obvious oscillatory behavior, the periods appear 
to be on the order of a few minutes. Thus, to improve the fits, I 
have accumulated the BCS data for longer intervals by requir- 
ing that the minimum number of counts in a channel exceed 
10,000 and that the accumulation time be at least 20 s. Note 
that since for most BCS data sets the minimum accumulation 
time is 3 s, this means that the data analyzed in this study are 
generally accumulated for at least 21s. The top panels in Fig- 
ure |3] show the Doppler shift data for the S XV and Ca XIX 
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channels for the flare data shown in Figure [2] accumulated in 
this manner along with the best-fit background trend. For both 
channels, the best-fit background was a two-term polynomial. 

The background-subtracted data along with the best-fit 
damped sine wave are shown in the bottom two panels of Fig- 
ure |3 In each wavelength channel, only the regions between 
the vertical dashed lines were fitted. While that data accu- 
mulated for longer times are still noisy, the fits appear to be 
reasonable. For the S XV and Ca XIX channels the reduced % 2 
values are 5.8 and 1.2, respectively. 

Table^summarizes the results of fitting all the flares exam- 
ined in this study for which reasonable fits to a damped sine 
wave were possible. For each flare, the table lists the date and 
time of the first S XV spectrum that was included in the fit, 
the class determined from GOES observations, the flare lo- 
cation on the Sun, and the fit parameters and the 1 a error in 
each parameter. Flare locations were taken from the NOAA 
Space Environment Center daily reports. If no position was 
available from that data source, the location of the flare was 
determined using Yohkoh Soft X-Ray Telescope (SXT) data. 
Those locations are marked with a 'y' in the table. For each 
flare, there are two rows in the table. The first row is for the 
S XV Doppler shift measurements and the second row is for 
the Ca XIX measurements. Although the origin of time for 
the fits in each channel was tied to the first fitted spectrum, 
the values for the Ca XIX fit parameters in the table have been 
adjusted to correspond to the initial time used for the corre- 
sponding S XV fit. This makes the fitting parameters for the 
two BCS channels directly comparable. The errors listed for 
each parameter in the table are the diagonal elements in the 
covariance matrix for the fits. 

Often there are significant Doppler-shift fluctuations before 
the time I have selected to begin fitting the data. The impres- 
sion one gets looking at the data is that the initial Doppler- 
shift disturbance is somewhat distorted — making it difficult 
to fit a decaying sine wave to both this initial fluctuation and 
the later more regular signal. 

Examination of the table entries for individual flares shows 
that there are events for which all the fit parameters for the 
two BCS wavelength channels agree to within the errors and 
events for which the fits result in one or more of the fit param- 
eters differing significantly. The fit parameters for the 1991 
October 21 event shown in Figures|2]and[3]are typical of many 
of the table entries. The frequencies and phases for the two 
channels are close in the two channels, while the amplitudes 
and decay rates are somewhat different, but still close enough 
to suggest that both channels are observing the same oscilla- 
tion. For some events, though, most of the fit parameters differ 
significantly, suggesting that the two BCS wavelength chan- 
nels are sampling different oscillations. The 1993 September 
26 flare is an example. 

Since the BCS images the entire Sun and the S XV and 
Ca XIX wavelength channels are sensitive to different temper- 
ature plasmas, it is certainly possible for each channel to be 
sampling different structures in the flaring region. One way to 
determine whether this is the case is to produce scatter plots 
of the fit parameter in the S XV channel against the same fit 
parameter in the Ca XIX channel for those flares where data 
from both channels could be fitted. Figure |4] shows those 
plots. With only a few exceptions, the periods, amplitudes, 
and phases, measured in the two BCS wavelength channels 
suggest that both channels are sampling plasma from the same 
structure. There is significantly more scatter in the decay time 
plots, but that parameter also has large errors associated with 



it. The errors are large for the decay time determination be- 
cause, as Figure [5] shows, generally only one cycle or less of 
the oscillation is observed. 

Figure [5] shows histograms of the distributions of the fit 
parameters and the inferred displac ement ampl i tude fo r the 
events listed in Tabled Following Wang et al. (2003a), we 
define the maximum displacement amplitude for each fit as 

A=A Q /(co 2 + X 2 )^ 2 . (3) 

The distributions shown in the figure appear to be similar for 
both BCS wavelength channels. Table|5]lists average values 
and standard deviations for the parameters in each channel 
and the combined data. Also listed in the table are the av- 
erage p roperties fo r oscil lations observed with TRACE from 
l Aschwanden et all (12002) and observed with SUMER from 
IWang et all (l2003al) . A reanalysis of the SUMER da ta gives 
similar values to those in the table (Wa ng et alJ20 05). 

4. BCS OSCILLATION MODES 

The transverse loop oscillations observed with TRACE 
are thought to be fas t kink mode MHD oscillations 
( lAschwanden et alJ I20021 Those observed with SUMER 
have been interpret ed as slow standing mode MHD waves 
(Wangetal. 2003b). Comparison of the averages in Ta- 
ble with the TRACE and SUMER results suggests that the 
BCS values are closer to those from TRACE than those from 
SUMER. The overall picture, however, is not entirely clear. 
Because exposure times for the SUMER observations tended 
to be quite long (> 50 s), those observations are not sensitive 
to the shortest periods observed with the BCS. On the other 
hand, as shown in Figure |2j the BCS data often only extend 
for 20 minutes or less. Thus, for many events, they would not 
be sensitive to the longer periods seen in the SUMER data. It 
might be speculated, though, that the slow Doppler shift evo- 
lution seen in the S XV data in Figure|2]and partially removed 
by the background polynomial fits to the data is actually a 
portion of a longer period damped oscillation that is superim- 
posed on those being measured in this paper. 

For a loop of length L, the period of a slow-mode wave is 
given by 

2L 

P=—, (4) 

JC T 

where j is the node number, usually assumed to be 1, and cj is 
the slow magne toacoustic speed, which is close to the sound 
speed dRoberts et alJ l984). Early in the time period where the 
BCS Doppler shift oscillations are observed, the temperatures 
measured in the S XV channel are around 12 MK, while those 
measured in the Ca XIX channel are around 14 MK (see, e.g., 
Figure|2}. For a coronal composition plasma, the sound speed 
is given by c s w 0.152\/7 7 km s _1 . Thus, if the oscillations 
are slow standing mode waves, the average periods in the table 
and their standard deviations imply loop lengths of about 82 ± 
35 Mm for S XV and 99 ± 55 Mm for Ca XIX— roughly 113 
and 137 arcsec, respectively. For semicircular loops with the 
above lengths, this corresponds to loop radii of 26 and 32 Mm, 
respectively, or 36 and 44 arcsec. 

The plasma observed with the BCS S XV and Ca XIX chan- 
nels tends to have higher temperatures than those measured 
for the sam e events using t he Yohkoh soft X-ray telescope 
(SXT) (e.g., lDoscheklfT99a) . Thus, it is not clear that the 
the sizes of the loops observed with SXT are appropriate for 
comparison with the above estimates. Moreover, many of the 
flares are at the limb and have portions of the loop occulted 
by the solar disk. 
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FIG . 3 . — Averaged BCS Doppler shift data and decaying sine wave fits for the S XV and Ca XIX observations for the 1 99 1 October 2 1 flare. The zero value 
for the shifts plotted in the top panels is based on an arbitrary reference system for the BCS in which a flare near the solar Equator with no Doppler shift would 
have a value of mA. Time is measured from 02:47:34 UT. 




FIG. 4. — Plots of the fitting parameters for events measured in the S XV wavelength channel versus the Ca XIX parameters for the same flare. Error bars for 
each point are the lo~ errors for the fit. If both channels are observing the same event, then the data should scatter around the solid lines. 



There are, however, some interesting events. The 1992 July 
17 flare was not occulted and clearly showed a loop-like struc- 
ture. This flare was included i n a study of occu lted and nonoc- 
culted flares by Mariska & McTiernan ( 1999). The loop is 
near the limb, so projection effects make a length determina- 



tion difficult. A. Winebarger (2005, private communication) 
has used the SXT Be filter data for this flare to examine pos- 
sible geometries for the loop. The geometry is difficult to 
define, but if the aspect ratio is near one (a semicircular loop), 
the length would be approximately 60 Mm. The frequencies 



6 



Mariska 



10 



a> 6 




15 



8 
6 

4 - 
2 

0_ 





Period (min) 



5 10 
Decay Time (min) 



15 



E 10 




Amplitude (km s 



80 2000 4000 6000 8000 10000 
Displacement (km) 



FIG. 5. — Histograms showing the distributions of period, decay time, amplitude, and displacement measured from the BCS data. Each bar shows the total 
number of events in the bin with the dark gray showing the S XV contribution and the light gray showing the Ca XIX contribution. 



listed for this flare in Table [pimply lengths of 47 and 61 Mm 
for the S XV and Ca XIX observations, respectively. Given the 
uncertainties in connecting the SXT images to the plasma be- 
ing observed with the BCS S XV and Ca XIX channels, these 
numbers are in excellent agreement. 

If the oscillations are slow standing mode waves, other ob- 
servable physical quantities should also be oscillating. For 
a sound wave, v = c s 8p/p. Thus, for a 12 MK plasma, 
dp/p w 0.04, a very small density fluctuation. Even though 
the intensity is proportional to the square of the density, this 
is a small fluctuation. Also, since the BCS images the en- 
tire loop, it would be difficult to see this small an oscillation 
when the entire loop is observed since for the fundamental 
mode the density oscillation has anti-nodes at the two foot- 
points of the loop, and the oscillations at these two anti-nodes 
are in anti-phase. For the partially occulted flares, one might 
expect, however, to see intensity fluctuations. 

The four panels of Figure |6] show the results of fitting 
a polynomial background and exponentially decaying sine 
wave to the intensity data for the 1991 October 21 flare. Only 
the data between the vertical dashed lines were fitted. The 
top panels show the intensity data and the best-fit background 
trend. In both BCS channels, the background was fit with a 
cubic polynomial. Error bars are not included in the plots. 
Since the integration times for the portion of the data used 
in the fits were generally 21s, the errors for a count rate of 
2500 counts s _1 are about 11 counts s . For a 1500 counts 
s~ 1 count rate, the corresponding error is about 8 counts s~ 1 . 
Thus the fluctuations fitted in the bottom two panels are well 
above the level expected from counting statistics. 

Table[5]summarizes the intensity fitting results for this flare. 
While the origin of the time values for the fitting calculation 



in each channel was tied to the time of the first fitted spectrum, 
the results in the table for the Ca XIX channel have been ad- 
justed to correspond to the initial time used for the S XV fits. 
The beginning time for the S XV intensity fits is the same as 
the beginning time used for the velocity fits for this flare listed 
in Tabled Thus the intensity results for the two channels can 
be directly compared with each other and with the results for 
the same flare shown in Tabled 

The periods and decay times of the intensity fluctuations are 
comparable with the characteristics of the fits to the Doppler 
shift measurements for this flare. The phases for the fits to the 
intensity measurements are, however, different. This is most 
easily visualized by comparing the times of the first peak in 
the fits to the Doppler shifts in Figure with those of the fits 
to the intensity fits in Figure [6] For S XV, the intensity peak 
occurs 1 .54 minutes after the Doppler shift peak. This corre- 
sponds to a delay of 1/4 of the period, strongly suggesting that 
the oscillations observed in this flare are due to slow standing- 
mode MHD waves JSakurai et alJ2002l) . 

Additional evidence for intensity oscillations in the 1991 
October 21 flare is provided by data from the GOES satellite. 
One-minute integration time data for this flare are available 
from the National Geophysical Data Center. Polynomial fits 
to the same time intervals used for the S XV intensity data for 
this flare reveal the same oscillations in the residual signal in 
both the 1-8 and 0.5-4 A GOES channels. 

One curious feature of the intensity fluctuations appears to 
be a tendency for the oscillation in the residual signal to be 
more obvious in the data from higher temperature plasma. 
Note, for example in Figure [6j that the oscillations in the 
Ca XIX intensity are better defined than those in the S XV 
residuals. This also appears to be the case for the two GOES 



Characteristics of Doppler Shift Oscillations 



7 



3000 
2500 
2000 
1500 









SXV 






- X 










5 10 
Time (min) 




5 10 15 20 
Time (min) 



FIG. 6. — Averaged BCS S XV and CaXIX intensity data and decaying sine wave fits for the 1991 October 21 flare. Only the data between the vertical dashed 
lines have been fitted. Time is measured from 02:47:34 UT — the same starting time used for Figurel3l 



channels. Examination of the the intensity data from the BCS 
Fe XXV channel also shows oscillations in the residual inten- 
sities for this flare that are well defined. This is probably due 
to the reduced amount of background emission at higher tem- 
peratures that must be removed from the oscillatory signal and 
suggests that oscillations should be detectable in the lower en- 
ergy channels for occulted flares observed with the Hard X- 
Ray Telescope on Yohkoh and those observed with the Ramaty 
High Energy Solar Spectroscopic Imager. 

Intensity fluctuations are not present in all the flares that 
exhibit Doppler shift oscillations. Thus it is worthwhile to 
ask whether some of the flares might be exhibiting fast kink 
mode oscillations which correspond to the displacement of an 
entire flux tube. For a loop of length L, the period of the wave 
is given by 



where j is the node number, usually assumed to be 1, and 
q is the kink speed (Roberts et al.l ll984l) . For a flux tube 
with an internal magnetic field strength equal to the external 
field strength, but a much greater plasma density, « V2ca, 
where ca is the Alfven speed. 

For the 1992 Jul 17 event discussed above, the estimated 
loop length of 60 Mm and a period of 3.3 min (the average of 
the S XV channel and Ca XIX channel results for this event) 
result in an Alfven speed of about 600 km s . Electron 
densities from 10 10 to 10 12 cm then imply magnetic field 
strengths ranging from 27 to 275 G. For a number of flares 
observed with the Naval Res earch Labo ratory S082A spec- 
trograph on Sky lab, Keenan et al. ( 1998) used density sensi- 
tive emission lines of Ca XVI to determine an average electron 
density of 2.9 x 10 11 cm~ 3 . This value would imply a mag- 
netic field strength of about 128 G. The low end of this range 
of field strengths is consistent with ranges that are typically 
quoted for flaring loops (e.g., Tsuneta 1996). 

Taking a temperature of 14 MK, a density of 3 x 10 1 1 cm -3 , 



and magnetic field strength of 128 G results in a plasma j8 of 
about 1.8. For a density of 10 10 cm~ 3 and a field strength of 
27 G, the plas ma /3 is about 1.3. The same conclusion was 
reached by Wang et al. (2002), who obtained values near two 
for j3 . Generally, the coronal plasma is thought to be low beta 
plasma. Moreover, equation Q has been derived by assuming 
that the plasma is low-/3 . Thus further analysis of the possibil- 
ity that some of the observed Doppler shift oscillations are due 
to kink modes will require a more careful examination of the 
governing equations. The existence of intensity fluctuations 
shifted by 1 /4 phase in some of the flares argues strongly that 
all of the flare Doppler shift oscillations observed with BCS 
are slow mode standing waves. 

5. DAMPING OF THE OSCILLATIONS 

Linear slow wave dissipation theory predicts that the decay 
time for an oscillation should scale as P , where P is the oscil- 
lation period ( P orter et alll 19941). More detailed calculations 
predict different slopes. Ofman & Wans (2002) have inves- 
tigated the damping of the slow-mode oscillations observed 
with SUMER using one-dimensional hydrodynamic simula- 
tions. They considered both damping by thermal conduction 
and compressive viscosity and found that for the SUMER ob- 
servations the damping was dominated by thermal conduc- 
tion. For the conditions in the loops observed with SUMER, 
they found that the scaling goes as P lA1 for a loop temper- 
ature of_6 ; 3_J^K_and as P 1 for a loop temperature of 6.8 
MK. lOfman & Wang l J2002I) attributed the smaller slope to 
the nonlinearity of the observed oscillations; they assumed 
A /c s = 0.18. 

Figure shows the decay time of the oscillations plotted 
against the period for all the events observed with the BCS. 
The line in the plot is the best fit to all the data and cor- 
responds to f(j eca y = 2.30P ' 39 — a much smaller slope than 
that determined by Wang et al. (2003a) from the SUMER ob- 
servations. The BCS results, however, apply to a much higher 
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FIG. 7. — The measured decay time vs. period for the events observed with 
the BCS. Error bars (Iff) for both the period measurement and the decay 
time measurement are shown. The solid line shows the best fit relationship 
between the log of the period and the log of the decay time. 



temperature plasma, 12 to 14 MK, than the SUMER observa- 
tions. (While the peak in the S XV emissivity function is at 
about 6.8 MK, the spectral fitting to the observations, which 
include temperature-sensitive lines, results in the higher tem- 
perature.) For the SUMER data, the highest temperature 
emission line in which oscillations were typically observed 
was from Fe XXI with a f ormation temperature of about 7.0 
MK. The Ofman & Wang| (12002ft calculations show that the 
slope decreases for increasing plasma temperature. Since the 
damping is by conduction, which scales as T 25 , one might 
expect a sizable change in the slope as the loop temperature 
increases. On the other hand, the oscillations observed with 
the BCS have a much smaller amplitude. For a temperature of 
12 to 14 MK, the average amplitude of 17.1 km s _1 leads to 
/ c s ~ 0.04, somewhat less nonlinear than is the case for the 
SUMER oscillations. Moreover, there is clearly considerable 
uncertainty in the BCS results, especially in the decay time 

determination; 

While the lOfman & Wangl (12002) calculations suggest that 
comparisons between the BCS and SUMER observations may 
be hindered by the different formation temperatures of the 
emission lines being observed, it is useful to plot all the 
data together, especially since the BCS observations explore 
shorter periods than are available with the SUMER data. Fig- 
ure|8]shows the decay time vs. period for the BCS events plot- 
ted in Figure0along with the results for events observed with 
SUMER and TRACE For the SUMER data, only the 49 cases 
fitted by Wang et al. (2003a) have been included. Note that 
while the TRACE data points lie in the same period interval on 
the figure as the BCS data, they generally have longer decay 
times than the BCS data — further suggesting that the TRACE 
data represent a different kind of oscillation. Of course, the 
loops measured in the TRACE observations contain plasma 
at significantly lower temperatures, ~ 1 MK, than those ob- 
served with the BCS and SUMER. 

Also plotted on the figure are the least squares fits to the 
BCS data that are shown on Figure0 the SUMER data, and 
the combined BCS and SUMER data. The fit to the combined 
BCS and SUMER data results corresponds to the expression 
= 1.30P - 81 . 



ecay 




6. DISCUSSION 



It is generally assumed that the emission observed with 
the BCS is produced in the loop or loops that have been 
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FIG. 8. — The measured decay time vs. period for t he even ts observed with 
the BCS. the SUME R events reported bvlWang et al J 12003d). and the TRACE 
events reported by Ofman & Aschwanden 1 2002 ), Wane & Solanki 1 2004 ), 
and Verwichte et al. i 2004|). The dashed line shows the best fit relationship 
between the log of the period and the log of the decay time plotted in FigurelTl 
the dot-dash line shows a fit to the Wang et al. 1 2003a) data, and the solid line 
shows the fit to the combined BCS and SUMER data sets. 



heated by the flare energy release process. Thus the ob- 
served oscillations are probably triggered by the heating due 
to the flare energy release. The observation of damped slow- 
mode oscillations in the high-temperature flare plasma then 
presents interesting challenges for some of the current mod- 
els for the evolution of flaring loops. TRACE observations 
clearly show that in the decay phase solar flares exhibit many 
loop-like structures (e.g.. [ Warren et alJll999l iWarrenlfeOOOt 
Aschwanden & Alexander 2001). This has led to attempts to 
model a flare as a succession o f independently heated threads 
(Warr en & Doschekl2005l:IWarrenll2005l) . To model a GOES 
M2.0 flare. lWarrenT(i2005l) needed 50 individual threads, with 
each one introduced 50 s after the previous thread. Each 
thread was heated using a triangular heating function with a 
width of 200 s. 

The model fits the BCS S XV and Ca XIX channels well, yet 
over the roughly 5-10 minute period where one sees oscilla- 
tions in most of the events discussed here, many individual 
threads contribute to the signal observed in the BCS chan- 
nels. This would seem to suggest that either multiple threads 
contain damped slow-mode waves with the same parameters, 
or that an individual thread is dominating the emission ob- 
served over the time interval whe re the o scillations are ob- 
served. The model described in Warren (2005) is still un- 
der development, so it is not clear, for example, what the 
time dependent Doppler shifts look like when the ensemble 
of threads is considered in their entirety as might be the case 
for a flare observed on the disk, or when only the top portions 
of the threads are considered as might be the case for a flare 
that is partially occulted by the solar limb. It is also not yet 
clear what the minimum number of threads is that will still 
smoothly reproduce the light curves observed with the instru- 
ments with which the data have been compared. 

Yohkoh observe d the flares that triggered a number of the 
events analyzed by Aschwand en et alJ (12002). Only the 1999 
October 25 flare exhibited oscillations in the BCS data that 
could be fit using equation Q. The three features observed 
with TRACE exhibited oscillation periods that ranged from 
2.38 to 2.70 minutes, with one feature showing a decay time 
of 3.33 minutes. In contrast, the BCS S XV and Ca XIX obser- 
vations show periods of 9.7 and 9.1 minutes and decay times 
of 10.0 and 5.68 minutes, respectively. Thus, while the os- 
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cillations may be triggered by the same event, it is unlikely 
that TRACE and the BCS are observing the same oscillations. 
This is not surprising given the significant differences in tem- 
perature sensitivity between the BCS and TRACE. 

Yohkoh SXT data are available for 18 of the 20 flares in- 
cluded in Tabled Most of the flares are at the limb and very 
little looplike structure ca n be discerned in th e data. Ten of the 
flares are included in the iSato et alJ JT991 catalog of flares 
observed with the Yohkoh SXT and Hard X-Ray Telescope 
(HXT). The images in the catalog made from the HXT data 
generally show only one emitting feature — again making it 
difficult to determine loop size parameters. In none of the 
flares have I been able to determine whether the oscillations 
are taking place in a pre-existing loop or the flaring loop it- 
self. Some of the SXT data sets, however, provide extensive, 
high-cadence observations of the event. Studies of the time 
histories of the intensities in the individual pixels may pro- 
vide additional insight into the nature of the oscillations. This 
will be the topic of a future investigation. 

Of the 103 flares observed by the Yohkoh BCS for this study, 
only the 20 in Table^exhibited Doppler-shift oscillations that 
could be easily fit with a well-defined decaying sine wave. As 
I mentioned earlier, an additional 18 events could be fitted, 
but were not used. Many of the other events, though, did ex- 
hibit significant Doppler-shift fluctuations. Since the BCS ob- 
serves all the flaring plasma at the temperature of formation 
of the S XV and Ca XIX lines, this is probably an indication 
that multiple structures are exhibiting Doppler shift fluctua- 
tions that are out of phase. While no attempt was made to 
observe flares at all center-to-limb positions on the disk, the 
results in Table [l] clearly suggest the flares near the limb tend 
to show well-defined oscillations more frequently than those 
closer to disk center. Since the Yohkoh BCS observes the en- 
tire Sun, this may simply imply that flares for which a portion 
of the flaring loops are occulted by the solar limb provide a 
cleaner Doppler-shift signal. 

Fitting the intensity fluctuations using equation Q is more 
difficult than fitting the Doppler-shift fluctuations. While the 
background trend in the Doppler shift data can usually be re- 
moved with just one or two terms in equation Q, the underly- 
ing trend in the intensities usually requires more terms. If one 
tries to detrend the intensity data with the peak in the inten- 
sity included, the background removal is even more complex. 



Moreover, including the background terms as adjustable fit- 
ting parameters often results in the fit selecting a damped sine 
wave with a substantial amplitude, up to 25% of the peak in- 
tensity, with a significantly altered background light curve, 
which is probably not a true representation of the oscillation. 

While I will defer a detailed analysis of the intensity fluctu- 
ations to a future paper, a preliminary examination of the S XV 
intensity fluctuations for the 17 flares in Table [J for which 
S XV Doppler-shift fits were available was conducted. All 17 
flares could be fitted with expressions of the form of equa- 
tion Q. Of those 17, only six showed phase lags that were 
judged to be consistent with that expected for a slow standing- 
mode MHD wave. Fits to five of the remaining flares resulted 
solutions that were viewed as having unrealistic background 
terms. Fits for the remaining six flares looked reasonable, but 
the phases were not close to the expected 1/4 period lag. Five 
of the six flares with reasonable phase lags were at the limb. 
The flares at the limb generally show only a small concen- 
trated emitting region in the SXT data. 

Detection of damped intensity fluctuations for some flares 
in the S XV, Ca XIX, and Fe XXV wavelength bands and in 
broad-band GOES observations raises the intriguing possibil- 
ity that spatially-resolved coronal loop observations with suf- 
ficient signal-to-noise and short enough exposure times may 
provide a wealth of information on the standing slow-mode 
waves. It may be possible to d etect dam ped int ensity oscil- 
lations in the Yohkoh SXT data. iMcKenzie & Mullanl J 19971) 
have searched for oscillations in SXT image sequences, find- 
ing some evidence for relatively short period fluctuations, 9.6 
to 61.6 s. Their time-series analysis may, however, have 
missed the kinds of intensity fluctuations reported here. The 
high sensitivity and spatial resolution Soft X-Ray Telescope 
and EUV Imaging Spectrometer instruments on the upcoming 
Solar-B satellite scheduled for launch in 2005 August should 
prove ideal for this kind of observation. 
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TABLE 1 

Doppler Shift Oscillation Properties 



to (UT) 




Class 


Location 


A 
(fans- 1 ) 


m 

(rad min~') 



(rad) 


X 

(min -1 ) 


1991 Oct 7 


10:16:57 


C9.9 


S13W82 


66.18±5.91 


1.36 ±0.04 


-5.80 ±0.07 


0.273 ±0.047 










38.92 ± 11.54 


1.27±0.10 


—5.63 ±0.06 


0.653 ±0.138 


1991 Oct 21 


02:47:34 


C4.3 


S06E89y 


18.89±1.04 


0.99 ±0.04 


-5.35±0.17 


0.276 ±0.020 










10.82 ±2.22 


1.17±0.12 


-5.54±0.17 


0.497±0.131 


1992 Jan 15 


08:59:21 


C5.0 


S18W90 


6.93 ±0.76 


1.39±0.09 


-5.54±0.17 


0.388 ±0.086 










9.69 ±2. 10 


1.10±0.08 


-4.98 ±0.30 


0.253 ±0.068 


1992 Feb 26 


01:35:57 


M1.3 


S15W90 


22.63 ±2.74 


2.47 ±0.11 


-5.82±0.21 


0.498 ±0.056 










20.40 ±6. 11 


1.93±0.10 


-4.00 ±0.27 


0.298±0.115 


1992 Mar 16 


06:06:15 


C5.8 


S25W21 


9.09 ±0.28 


1.32 ±0.05 


-0.84 ±0.04 


0.404 ±0.028 










7.93±2.17 


1.44 ±0.07 


-0.24 ±0.20 


0.160±0.110 


1992 Apr 19 


02:11:32 


C3.9 


N15E78 


6.49 ±1.74 


2.16±0.05 


-0.50 ±0.20 


0.1 17 ±0.063 










4.38 ±2.96 


1.58±0.23 


-5.17±0.59 


0.193 ±0.285 


1992 Apr 20 


06:33:25 


C3.6 


S15E39 


5.89±1.05 


1.94 ±0.05 


-4.69 ±0.21 


0.149 ±0.043 


1992 May 27 


14:37:53 


C7.1 


N23W59 


11.01 ±2.05 


2.48 ±0.06 


-3.58±0.13 


0.391 ±0.090 








6.00 ±1.84 


2.22 ±0.07 


-2.30±0.26 


0.1 17 ±0.075 


1992 Jul 4 


22:49:16 


C6.9 


S12E89 


6.22±1.93 


0.67 ±0.06 


-2.94 ±0.13 


0.229 ±0.078 










6.45± 1.10 


0.58 ±0.08 


-2.74 ±0.22 


0.181 ±0.098 


1992 Jul 17 


22:36:31 


C5.3 


S11W87 


4.49 ±1.23 


2.09 ±0.05 


-2.62 ±0.17 


0.214 ±0.079 










7.87±4.82 


1.76±0.11 


-1.16±0.52 


0.202±0.139 


1992 Aug 24 


01:12:48 


C2.2 


N14W90 
















7.58 ±2.94 


2.99 ±0.25 


-2.90 ±0.40 


0.474 ±0.217 


1992 Oct 12 


21:50:44 


C2.5 


S16W88 


11.37±1.87 


1.60±0.06 


-2.43 ±0.05 


0.284 ±0.090 










8.77 ±2.79 


1.76±0.17 


-2.86 ±0.41 


0.235 ±0.130 


1992 Oct 27 


22:14:55 


C5.4 


N08W90 


14.03 ±1.44 


0.96 ±0.04 


-1.07 ±0.20 


0.165 ±0.028 










16.96±3.94 


1.25 ±0.14 


-2.18±0.46 


0.156 ±0.107 


1992 Nov 5 


06:24:08 


M2.0 


S17W88 


9.22 ±0.60 


0.92 ±0.02 


-4.37 ±0.11 


0.072 ±0.017 










24.41 ±2.00 


0.74 ±0.02 


-2.90±0.10 


0.159±0.018 


1993 Feb 1 


10:07:38 


C6.9 


S11E90 


3.65 ±0.32 


1.55 ±0.23 


-1.94 ±0.54 


0.270 ±0.058 










3.52±1.66 


1.69 ±0.26 


-1.64 ±0.55 


0.253 ±0.219 


1993 Apr 15 


09:08:36 


C1.2 


S19W90y 


26.00 ±1.52 


0.78 ±0.01 


-1.60 ±0.05 


0.142 ±0.009 






42.19±4.08 


0.83 ±0.03 


-2.04 ±0.14 


0.193 ±0.016 


1993 Sep 26 


10:22:47 


C3.4 


N12E90 


15.43 ±4.23 


0.96 ±0.03 


-3.69 ±0.25 


0.173 ±0.030 








76.32± 17.88 


0.42 ±0.03 


-0.68 ±0.04 


0.193 ±0.053 


1994 Feb 27 


09:11:09 


M2.8 


N08W88 


3.60 ±0.60 


1.39±0.03 


-0.55 ±0.04 


0.143 ±0.031 










4.75 ±1.16 


1.40 ±0.04 


-0.67±0.19 


0.190 ±0.055 


1999 Oct 25 


06:36:18 


M1.7 


S18E88y 


17.25 ±1.52 


0.65 ±0.02 


-0.10±0.09 


0.100±0.017 








33.17±4.29 


0.69 ±0.02 


-0.35 ±0.06 


0.178±0.037 


2000 Jul 14 


00:41:45 


M1.5 


N20W80y 


















43.01 ±7.50 


0.96 ±0.04 


-5.88±0.12 


0.341 ±0.057 



TABLE 2 

Average Doppler Shift Oscillation Properties 



Parameter 


SXV 


CaXIX 


Combined 


TRACE 


SUMER 


Oscillation Period (min) 


5.2±2.2 


5.8 ±3.2 


5.5±2.7 


5.4±2.3 


17.6±5.4 


Decay Time (min) 


5.5±3.1 


4.7 ±1.8 


5.0±2.5 


9.7 ±6.4 


14.6±7.0 


Amplitude (km s -1 ) 


14.4±14.5 


19.6± 19.2 


17.1 ±17.0 


42 ±53 


98 ±75 


Displacement (km) 


730 ±742 


1400 ±2260 


1070±1710 


2200 ±2800 


12500 ±9900 


Decay Time to Period Ratio 


1.13 ±0.65 


0.98 ±0.62 


1.05 ±0.63 


1.8 


0.85 ±0.35 



TABLE 3 
Intensity Fitting Results 



Parameter S XV Ca XIX 



A (counts s- 1 ) 41.5±8.0 38.8±7.1 

o (rad min"') 1.05 ±0.04 0.81 ±0.04 

<|>(rad) -0.72 ±0.25 0.53 ±0.25 

A (min- 1 ) 0.12± 0.04 0.14 ±0.02 



